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A recent paper reported elliptically polarized high-order harmonics from aligned N2 using a lin- 
early polarized driving field [X. Zhou et al, Phys. Rev. Lett. 102, 073902 (2009)]. This observation 
cannot be explained in the standard treatment of the Lewenstein model and has been ascribed to 
many-electron effects or the infiuence of the Coulomb force on the continuum electron. We show 
that non-vanishing ellipticity naturally appears within the Lewenstein model when using a multi- 
center stationary phase method for treating the dynamics of the continuum electron. The reason 
for this is the appearance of additional contributions, that can be interpreted as quantum orbits 
in which the active electron is ionized at one atomic center within the molecule and recombines 
at another. The associated exchange harmonics are responsible for the non-vanishing ellipticity 
and result from a correlation between the ionization site and the recombination site in high-order 
harmonic generation. 

PACS numbers: 42.65.Ky, 42.65.Re 



I. INTRODUCTION 

High-order harmonic generation (HHG) is a highly 
non-Unear process in which a medium in an intense laser 
field emits coherent radiation at multiples of the driving 
frequency. The short time scale of the generation process 
makes HHG a promising source of coherent attosecond 
pulses in the XUV regime. When the target is isotropic, 
such as a gas of atoms or unaligned molecules, it follows 
from symmetry arguments that the emitted harmonics 
have to be polarized parallel to the polarization axis of a 
linearly polarized driving laser. Breaking this isotropy by 
using a sample of aligned molecules for HHG allows for 
a non- vanishing perpendicular polarization component. 
This scenario is sketched in Fig. [TJ The perpendicular 
component is generally heavily suppressed compared to 
the parallel component. This has lead to a focus in the lit- 
erature on the parallel component. Recent experiments 
on aligned N2, O2 and CO2, however, have reported a 
non- vanishing perpendicular component [ij. Moreover, 
elliptically polarized harmonics have been measured from 
aligned N2 and CO2 01 • The presence of elliptically po- 
larized harmonics opens up the possibility of generat- 
ing elliptically polarized attosecond pulses in the XUV 
regime. 

From a theoretical point of view, the observation of el- 
liptically polarized harmonics is very interesting because 
it serves as an important benchmark for different mod- 
els. Several approaches are currently used to calculate 
the HHG response. Ideally, one should propagate the 
TDSE [3, 4]. This approach has been used for atoms 
and small molecules. An example is the prediction of 
comparable polarization components from aligned HJ 
near minima in the spectrum 0. For systems beyond 
HJ and H2 one commonly uses the semi-classical simple- 
man's model or the quantum mechanical Lewenstein 
model 01 to describe HHG. Both models provide a three- 



step picture of the process in which a single active elec- 
tron is brought into the continuum due to its interaction 
with the laser field, propagates in the field, and subse- 
quently recombines and emits radiation. 

The calculation of HHG spectra using the Lewenstein 
model is generally combined with a stationary-phase 
method to reduce the computational cost. As discussed 
in Sec. Ill B[ this procedure causes the resulting harmonics 
to be linearly polarized parallel to the polarization of the 
driving laser. This shortcoming has been attributed to 
the Lewenstein model itself 1,^. However, the standard 
stationary-phase method is not in general valid for molec- 
ular targets as it fails to take into account the spatial ex- 
tent of the molecule Hj. Instead, an extended stationary- 
phase method can be applied to evaluate the Lewenstein 
model, which leads to so-called exchange harmonics 
These are contributions to the total harmonic emission 
that are caused by quantum orbits describing an ioniza- 
tion event at one atomic center followed by recombination 
at another atomic center [see Fig. [2]. 

We show that the Lewenstein model does lead to el- 
liptically polarized high-order harmonics when exchange 
harmonics are included. Exchange harmonics rely on 
a correlation between the ionization and recombination 
events in HHG. This is different to the perpendicular har- 
monic component reported in [9| , which is obtained using 
scattering states for the continuum dynamics, thus caus- 
ing a spread in momentum of the returning electron. The 
derivation in 9, 10] uses the standard stationary-phase 
method in order to separate the three steps of the HHG 
process. Consequently, any correlation, in the sense of 
exchange harmonics, between ionization and recombina- 
tion sites is lost. Detailed comparison with experiments 
like that in [i| will serve to determine the importance of 
this correlation. 

The present paper is organized as follows. We briefly 
describe the Lewenstein model for aligned molecules in 
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We define the pulse by the vector potential 

Fq 

A{t) = — f{t)cos{LJot)ex, 

UJQ 



(1) 



where Fq is the maximal field strength, luq the angular 
frequency, f{t) the pulse envelope starting at t ~ and 
ending at i = T (see Sec. IIIip . and the polarization 
direction. The electric field is obtained from Eq. ([T]) as 
F(t) = -dtA{t). 



A. The Lewenstein model for aligned molecules 
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FIG. 1: (Color online) (a) The driving pulse propagates along 
the z-axis and is linearly polarized along the x-axis. The 
molecule is aligned at an angle of (p to the a;-axis in the xy- 
plane. The alignment at a given angle with respect to the 
external field results in harmonic emission with components 
Ex and Ey. The hatched structure indicates the HOMO 
(cf. Sec. Ill), (b) The polarization ellipse for a single har- 
monic. The angle rj of the major axis and the length b of the 
minor axis have been exaggerated for clarity. 



Sec. Ill Al The stationary-phase method is introduced 
in Sec. Ill Bl and the resulting harmonics shown to be 
linearly polarized parallel to the linear polarization of 
the driving laser. In Sec. Ill CI we present the extended 
stationary-phase method, and Sec. IIIII contains results 
demonstrating that such a treatment of the Lewen- 
stein model leads to elliptically polarized harmonics from 
aligned N2. Wc summarize our findings in Sec. IIVI 
Atomic units [h = e = trie = ao = 1] are used throughout 
unless stated otherwise. 



II. THEORY 

We treat the response of a single molecule to the elec- 
tric field F(i) of an intense femtosecond laser pulse. The 
laboratory frame is chosen such that the pulse propagates 
along the z-axis and is linearly polarized along the x-axis. 



The calculated spectrum Sn{uj) of the harmonic com- 
ponent along n depends on the choice of gauge and form. 
Following Tl], we use the length gauge for the inter- 
action Hamiltonian, and the velocity form to determine 
the dipole acceleration. Thus, given the dipole velocity 
Vdip(i), the signal is 



(2) 



When evaluating Eq. ^ it is necessary to include con- 
tributions from all occupied molecular orbitals A fT5|, so 



(vdip(i)) = E^^^^w)' 



(3) 



with weights Px fulfilling the normalization condition 
Fx = 1. Introducing Euler angles TZ = {(f), 9, x) for 
the molecular orientation in the laboratory frame, the 
alignment distribution G{TV) for a given molecule and 
alignment pulse can be calculated as in [l^. Follow- 
ing this results in a coherent average over fixed 
alignment angles 



(vA(t)) - / d7^G(7^)(vA(7^,^)). 



(4) 



The dipole velocity of a single perfectly aligned molec- 
ular orbital is calculated using the Lewenstein formula 



(VA(7^,^)) 



^g-^S.(k,M-r)p(^ - r) • dio„,A(7e, k, t - r) 



-|-complex conjugate. 



(5) 



where the nuclei are taken to be fixed during the driving 
pulse. This latter approximation is valid for heavy atoms 
such as nitrogen [H, [l^ . 

Equation ([5]) depends on the Fourier transform of the 
electronic wavefunction ip\ in both the ionization term 

dion,A(k, t-T)= ^Vk+A(t-r)V^A[7^, k + A{t - t)], (6) 

and the recombination term 



(7) 
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The vertical ionization potential \E\\ enters into the 
phase factor Sx(k, t,t — t) as 



The Fourier transforms in Eqs. ©-(IT]) then separate into 
terms each relating to just one atomic center. This im- 
plies that 



S.iK t,t^r)^ i [k + Ait"f dt" - E,r. (8) ^ ^ ^ ^^^^^^ ^ _ 



(13) 



B. Standard stationary-phase method 

The momentum integral in Eq. ^ describes the prop- 
agation of the active electron in the continuum. This 
integral is usually calculated using the stationary phase 
method, approximating it by a normalization factor times 
the value of the integrand at the stationary point kg of 
the rapidly varying phase. 

Identifying the total phase with S\{\i, t,t — T) gives the 
stationary-phase condition 

Vk5A(k,i,t-T)|^^j^^ =0, (9) 

which leads directly to the stationary momentum 



and 



1 



A{t")dt". 



(10) 



The standard stationary-phase method thus yields 



(vA(7^,^)) 



^e~'SM.t-r)^^^ - r) • dion,A(7e, k„ t - r) 
-l-complex conjugate. (11) 

Calculated spectra depend only weakly on the chosen 
value of e, which is set to 1 in the following. 

Equations pO)) - (fTT|) select the dominant contribution 
to the full momentum integral in Eq. It is conven- 
tional to interpret the result of such a stationary phase 
analysis in terms of quantum orbits (TtI . [isj . In this lan- 
guage, the active electron tunnels out, propagates paral- 
lel to the polarization of the laser, and recombines. 

The vector character of PT|) is determined by 
Vroc,A('^, kg, t). Referring to Eqs. ([7]) and pil)l . it is clear 
that the emitted harmonics are polarized parallel to the 
linear polarization of the driving laser. The Lewenstein 
model with the standard stationary phase method cannot 
give rise to ellipticities different from zero. 



C. Extended stationary-phase method 



The treatment in Sec. Ill Bl is oversimplified in the 



molecular case 



The reason for this is that di, 



and 



pick up phase factors related to the internu- 



clear separation. This is readily seen when the electronic 
wavefunction is expanded around the atomic centers as 



V'A(r) 



(12) 



V,.oc,A(k,t) 



/ J reCjA 



(k,i). 



(14) 



Each term in Eq. (|T3|) -(fT4 |) acquires an associated phase 
factor e»[k+A(t-r)]-R„, e*[''+^(*)]-^"/ , which has to be 
taken into account when performing the stationary-phase 
analysis. Here ri^ refers to the atomic center involved in 
evaluating di^^^ ^, and n/ refers to v"/^ ^. 

The above discussion means that Eq. ([5]) can be split 
into 



(vA(7e,t)) 



(v-"^(7e,t)) 



(15) 



The stationary-phase condition for ^v"'"^ (7?., t)) is 



/ (k + A{t")) dt" ~ (R„, - R„, 

Jt-T 

with associated stationary momenta 

kr'"^ = -- r A(t")di" - (R„, 



Rri 



0, (16) 



(17) 



Applying this multi-center stationary phase method to 
the integral in Eq. (O finally yields 



(VA(7^,^)) 



dT 



271 



3/2 



v:7,(7^,k-'^t) 



xe 



iS^iK' ' ■t,t-T) 



xF(t-r).d-^,(7e,kr 
+coniplex conjugate. 



(18) 



The terms in Eq. p8p can be interpreted as semi- 
classical orbits of electrons that are ionized at R„; and 
recombine at R„^ . The effect of using the multi-center 
stationary phase method is to properly include orbits 
from one atom to another. The resulting harmonics 
are referred to as exchange harmonics Q. A similar 
idea was explored in an earlier paper [Toj . and its above 
threshold ionization analogue is well-known in the liter- 
ature 0,[2l|,[22. 

The vector character of the emitted harmonics is de- 
termined by v"j"j^ in Eq. (|18p . In analogy to the 

standard stationary-phase method, v"j"j^ is parallel to 



As seen from Eq. p?|) . this means that 



exchange harmonics will contribute with a non- vanishing 
perpendicular component if the molecular axis has a non- 
vanishing projection orthogonal to the driving laser. The 
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FIG. 2: (Color online) Two orbits with ionization at the peak 
of the field and recombination two thirds of an optical cycle 
later. The atomic centers of the molecule are sketched as 
black dots to the right. The laser polarization is taken to 
be in the i-direction. The dashed (blue) line represents an 
orbit that gives rise to direct harmonics, i.e., ionization and 
recombination at the same center. A transverse component 
has been added manually to distinguish the return path. The 
dash-dotted (red) line is an orbit that gives rise to exchange 
harmonics, i.e., ionization at one center and recombination at 
another center. Note that the units on the two axis differ. 



FIG. 3: (Color online) Orientation dependence of harmonic 
17 from the Sag HOMO of N2 with equilibrium internuclear 
distance Ro = 1.0977 A. We use an 800 nm, 6 x lO" W/cm^ 
driving field. The envelope (for the vector potential) is trape- 
zoidal with three optical cycles turn-on and turn-off and five 
cycles constant amplitude, (a) Harmonic phase (b) Harmonic 
intensity (c) Ellipticity (d) Angle of polarization ellipse major 
axis with respect to molecular axis. In (a) and (b) the solid 
(blue) curves refer to the parallel polarization component and 
the dash-dotted (red) curves to the perpendicular component. 



conventional stationary phase method only includes the 
direct harmonics, which is why the standard treatment 
of the Lewenstein model fails to reproduce, even qualita- 
tively, the findings of P, . 

An example of the two types of orbits is presented in 
Fig. [21 The electron is ionized at one atomic center at 
the peak of the field and recombines two thirds of an op- 
tical cycle later. Laser parameters are the same as those 
used in Sec. IIIII below. The orbit that gives rise to direct 
harmonics has been plotted with a small artificial trans- 
verse displacement in order to separate the two parts of 
the electron motion. Similar orbits exist which describe 
ionization at the other atomic center. 

We note in passing that the full Lewenstein expression 
in Eq. ([5]) is independent of whether a single-center or 
a multi-center expansion is adopted. The distinction is 
purely an artifact of using a stationary-phase method to 
approximate the momentum integral. In this setting, a 
multi-center expansion is preferable as it easily allows the 
geometric phase factors to be extracted. 



III. RESULTS AND DISCUSSION 

We present here a proof-of-principle calculation on N2 
in order to illustrate that Eq. (fT5)) indeed leads to ellip- 
tically polarized harmonics when using a linearly polar- 
ized driving pulse. We assume only contributions from 
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FIG. 4: (Color online) Same as Fig. |3] but for harmonic 27. 



the Sag orbital, which is the highest occupied molecular 
orbital (HOMO). The HOMO is obtained using the stan- 
dard quantum chemistry software package GAMESS- 
US [Is] with a triple zeta valence basis set and diffuse 
i-shells. We assume perfect alignment, and perform no 
focal volume averaging. The alignment geometry and 
the polarization ellipse are shown in Fig. [TJ Based on 
a wavelet analysis [2J|, the r-integral in Eq. (fT5)l is re- 
stricted to allow only short trajectories. This means that 
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FIG. 5: (Color online) Absolute value of the harmonic ellip- 
ticity e as a function of harmonic order for selected alignment 
angles of N2. The nuclei are fixed at their equilibrium po- 
sition. The solid (blue) curve refers to </> = 40°, the dashed 
(green) curve to <^ = 50° , and the dash-dotted (red) curve to 
(j> = 60°. See the caption of Fig. [3] for laser parameters. 



the upper bound is set to 0.67 times an optical cycle. The 
driving field is taken to be an 800 nm pulse with 11 cy- 
cles in a trapezoidal envelope with 3 optical cycles for 
the linear ramp-up and ramp-down. The peak intensity 
is 6 X 10^^ W/cm2. 

We follow the polarization conventions in [2^. This 
means that the harmonic of angular frequency co is as- 
sumed to be a perfect plane wave 

(t) — fli cos(wt + (5i)e^ + a2 cos(cjt + (52)ei, (19) 
= o cos{ujt + (5o)ej;' ±6 sm{ujt + 6o)eyi , (20) 

where e^;, ey, ai, 02, e^'^eyi, a, b are defined in Fig. [T] (b). 
Phases Si {i — 0, 1, 2) are found by assuming that the har- 
monic phase is equal to that of the dipole velocity [2^ . 

The two physically relevant quantities are the elliptic- 
ity e, and the angle 77 that the polarization ellipse major 
axis makes to the polarization axis of the driving pulse. 
These are determined by 

e = ±- , < tan^Ve < - 21 



and 



tan(2?7) = . cos(J2 - <5i ) ■ (22) 

1 - (01/02) 



A. Equilibrium nuclear spacing 

In this subsection we will restrict ourselves to the equi- 
librium nuclear spacing of Rq = 1.0977 A. Fi gures[31[a)~ 
(b) and |lla)-(b) show data for harmonics 17 and 27, 



FIG. 6: (Color online) Same as Fig. O but with the nuclei 
fixed at 2Ro. 
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FIG. 7: (Color online) Same as Fig. (3] but for harmonic 57 
and with R = 2Ro . 



with solid (blue) curves referring to the component par- 
allel with the linear polarization of the driving field and 
dashed (red) curves to the perpendicular component. 
Panels (a) illustrate how the parallel polarization com- 
ponent of the harmonic emission changes its phase much 
faster than the perpendicular component. Comparison 
with panels (b) shows that the rapid change in phase is 
accompanied by a drop in harmonic intensity, causing an 
increase in the ratio of perpendicular to parallel harmonic 
intensity. 

We determine the ellipticity e and the angle t] of the 
major axis using Eqs. (|2T|) - (|22)) . Results are shown in 
panels (c) and (d). Both are seen to be very angle- 
dependent, taking non- vanishing values where the par- 
allel component has a rapid change in phase. This varia- 
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tion in phase changes the hehcity of the harmonics over 
a very short angle interval, while the angle of the major 
axis has a definite sign for a given harmonic. 

Figure [5] shows the absolute value of the ellipticity as a 
function of harmonic order for selected alignment angles. 
The angular dependence can largely be understood by 
comparing panel (b) in Fig. [3] with that in Fig. 01 The dip 
in the intensity of the parallel component is fairly narrow 
and moves to larger angles for increasing harmonic order 
in agreement with the two-slit interference formula [27| . 
The resulting ellipticity does not compare well with [3]. 
The experimental ellipticities are an order of magnitude 
higher and hardly have any angle-dependence. 

The presented calculation shows that elliptically po- 
larized high-order harmonics are predicted within the 
Lewenstein model. Improvements have to be made in 
order to compare the model directly with experimental 
data. One of these is to include the lower-lying molec- 
ular orbitals (HOMO-1, HOMO-2, etc.). The inclusion 
of these is expected to cause a dynamical minimum in 
the parallel component @, which would change the de- 
tailed structure. Taking into account the distribution 
over alignment angles G{TZ) might also prove to be im- 
portant. 

B. Extended nuclear spacing 

As discussed above, exchange orbits give rise to ellipti- 
cally polarized harmonics when the parallel polarization 
component is suppressed by the two-slit interference con- 
dition. Changing the internuclear separation will have an 
effect on this interference. Figure [6] shows results for N2 
with R = 2Ro. The plotted ellipticity is seen to peak 
at higher harmonic orders than was the case at R = Rq. 
There is no significant change in the magnitude of the ob- 
tained ellipticities. Figure [7] shows the detailed behaviour 
of the 55th harmonic. A second interference minimum 
appears at large alignment angles, which is consistent 



with the second order minimum in the two-slit model. 



IV. CONCLUSION 

Motivated by recent experiments we have ad- 

dressed the issue of understanding how elliptically polar- 
ized harmonics arise in HHG when using a linearly po- 
larized driving laser. It was shown that the Lewenstein 
model does not allow the emission of harmonics with a 
non-vanishing perpendicular component if the standard 
stationary-phase method is used to evaluate the contin- 
uum dynamics. A multi-center stationary phase method 
was adopted in order to take into account the spatial ex- 
tent of the molecule. This leads to additional terms in the 
Lewenstein model compared to the standard stationary- 
phase method. These were interpreted as quantum or- 
bits describing an ionization event at one atomic center 
followed by propagation in the continuum and recombi- 
nation at a different atomic center [1]. We illustrated 
that the total harmonic emission from N2 calculated us- 
ing this method is elliptically polarized near minima in 
the spectrum. 

The electron continuum was described using only 
Volkov waves, as opposed to the use of scattering 
states [21, indicating that the appearance of elliptically 
polarized harmonics is in part the result of having an 
extended target. Further work is required to ascertain 
the importance of the Coulomb interaction between the 
active electron and the molecular ion. 
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